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Understanding the nature of multipartite entanglement is a central mission of quantum informa- 
tion theory. To this end, we investigate the question of tripartite entanglement convertibility. We 
find that there exists no easy criterion to determine whether a general tripartite transformation can 
be performed with a nonzero success probability and in fact, the problem is NP-hard. Our results are 
based on the connections between multipartite entanglement and tensor rank (also called Schmidt 
rank), a key concept in algebraic complexity theory. Not only does this relationship allow us to 
characterize the general difficulty in determining possible entanglement transformations, but it also 
enables us to observe the previously overlooked fact that the Schmidt rank is not an additive entan- 
glement measure. As a result, we improve some best known transformation rates between specific 
tripartite entangled states. In addition, we find obtaining the most efficient algorithm for matrix 
multiplication to be precisely equivalent to determining the optimal rate of conversion between the 
Greenberger-Horne-Zeilinger state and a triangular distribution of three Einstein-Podolsky-Rosen 
states. 

PACS numbers: 03.67.Mn, 03.65.Ud 



One of the greatest discoveries in quantum physics [l| is 
that a multipartite quantum system can be in a so-called 
entangled state. There are an uncountable number of 
entangled states realizable by any quantum system and 
a natural question is how they are related to each other 
— specifically, if two given states can be converted to 
each other through local operations and classical commu- 
nications (LOCC), i.e. a protocol in which no quantum 
information is exchanged among the subsystems. Ideally, 
one would like to have an efficient way to decide whether 
one entangled state can be transformed into another via 
LOCC as this question has utmost physical relevance. By 
nature, entangled states are fragile and highly suscepti- 
ble to decoherence, a process in which quantum powers 
are lost. LOCC protocols describe the cheapest and ex- 
perimentally easiest ways to convert entanglement while 
minimizing decoherence among parties separated by ar- 
bitrary distances. 

As the notion of probability is inherent to quantum 
mechanics, the more natural question is with what prob- 
ability p can be converted into \tp) under LOCC? For 
p = 1, the LOCC transformation is called determinis- 
tic and for a general nonzero p, the protocol is called 
stochastic (SLOCC). Transformations of the latter form 

are written as \cj)) ^^'^^^^^ |^^^ Yot: bipartite systems, the 
problem is completely solved. Nielsen has provided nec- 
essary and sufficient conditions for whether two states are 
deterministically convertible Probabilistically, any 
bipartite state can be transformed into if and 
only if the matrix rank of the reduced density operator 
of |(/)) is greater than that of Furthermore, Vidal 



Q has derived a simple formula that gives the optimal 
probability for conversion. 

When the number of subsystems is greater than two, 
the situation becomes much more complicated. No longer 
can SLOCC convertibility be determined by examining 
the ranks of the reduced density matrices of the ini- 
tial and final states. For example, a system of three 
qubits can be partitioned into six equivalence classes 
defined by SLOCC convertibility between states in the 
same class However, two of the classes are indis- 
tinguishable by examining the ranks of each subsystem's 
density matrix. The states \GHZ) = ^(|000) + |111)) 

and \W) = ^(|001) + |010) + |100)) are representatives 
of each class respectively. Progress toward understand- 
ing SLOCC convertibility among four qubit states has 
also been made as the state space for these systems can 
be partitioned into nine different families of equivalence 
classes [Ij. Nevertheless, for both three and four qubit 
systems, the separation into SLOCC equivalence classes 
does not provide a solution for determining whether two 
states are SLOCC related simply because one must first 
determine to which classes the states belong. 

In this letter, we ask whether there is some relatively 
simple criterion for determining the convertibility of arbi- 
trary tripartite states like there is for bipartite states. As 
a complete solution to the convertibility problem should 
be able to determine whether one state can be trans- 
formed into another with a nonzero probability, we focus 
our attention on the class of SLOCC protocols to judge 
the difficulty of the complete problem. Ultimately we find 
that no simple criterion exists for testing the possibility of 
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a general tripartite entanglement transformation. In ad- 
dition, through the course of investigating this problem 
many other interesting results are obtained concerning 
specific tripartite transformation rates. The novel con- 
version rates are derived in part from our observation 
that the Schmidt measure (to be defined below) is not 
an additive quantity, something previously thought to be 
true 0. We now summarize our main findings: 

Denote by \^^) the unnormalized tripartite state where 
any two parties share an (unnormalized) EPR state |<I>) = 
|00)-t-|ll): 

1$'^) = \^)AB\'i>)Ac\^)BC 

= (|00).4|00)B + |10)A|10)i3)|00)c 



|00)A|01)B + |10)A|ll)i3j|01)c 

|01)a|00)s + |11)a|10)b)|10)c 
|01)a|01)b + 111)^111)5) |ll)c. 



(1) 



Theorem 1. 



(a) For any general tripartite and \ip), determining 
whether can be obtained from j?/;) by SLOCC is 
NP-hard. 

(b) \GHZ)®^ SLOCC ^^^^2^ 



(c) \GHZ)' 



SLOCC 
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(d) Let A = inf{u : \GHZ)®^^''-^ SLOCC |^3^^„ 
sufficiently large n}. Then A is precisely the ex- 
ponent for matrix multiplication, i.e., the smallest 
real number to such that two N by N matrices can 
be multiplied with 0{N'^) number of multiplications 
between linear functions on entries of the first ma- 
trix and linear functions on entries of the second 
matrix. 

Previously, only one copy of the W state is known to 
be convertible from three copies of GHZ and result (b) 
provides an improvement to this rate. Transformation 
(c) is important because it reveals that the three-party 
EPR extraction rate from GHZ is greater than one, a 
previously unknown possibility. Result (d) shows that 
existing lower and upper bounds for matrix multiplica- 
tion translate to lower and upper bounds on the optimal 
conversion rate between groups of the given states. This 
connection implies that finding the GHZ to three-party 
EPR conversion rate is highly difficult since the complex- 
ity of multiplying matrices is one of the most challenging 
open problems in computation theory. 

Our main technical tool is tensor rank, a key con- 
cept in algebraic complexity theory Q that has also 
been used to measure multipartite entanglement under 
the synonymous names of Schmidt rank and Schmidt 



measure [6[. The tensor rank of a multipartite state 
10) G Hi ® H2 'S) ■ ■ ■ (E) H„, denoted by rk{\(j))), is the 
minimum number r such that there exists G Hi, 

1 < j < r and 

r n 

j=i 1=1 

The quantity log2{rk{\(j))) is called the Schmidt measure 
of 10), denoted by sc/i(|0)). 

Tensor rank has been used in algebraic complexity the- 
ory as it captures the complexity of computing a set of 
bilinear maps and in particular the multiplicative com- 
plexity of multiplying two matrices. A set of bilinear 
maps are polynomials with respect to two distinct groups 
of indeterminates. The multiplicative complexity of the 
set is the minimum number of multiplications between 
the two groups required to evaluate all the polynomials. 
The multiplication of two N x N matrices produces a 
set of N'^ bilinear maps, one for each entry in the N x N 
product. The complexity of N x N matrix multiplication 
is denoted by fJ.{N, N) and the current best upper and 
lower bounds for fiiN,N) are 0{N^-^'^) and fiV^ - 3A^ 
respectively The complexity of matrix multipli- 

cation is also expressed as fi{N,N) = 0{N'^) where uj 
is called the exponent for matrix multiplication and de- 
fined as the smallest real number such that an algorithm 
exists for multiplying two N x N matrices using 0{N'^) 
multiplications. While uj is hypothesized to be two, de- 
termining the validity of this conjecture is a major open 
problem in computational science. For more details, a 
good reference is chapter 28 of [l3l- 

Tensor rank analysis has already shown to be valuable 
in quantum information as it is the distinguishing prop- 
erty between the \G HZ) and \W) equivalence classes of 
three qubits [1, [Tlj. It has also been useful in charac- 



terizing the entanglement in graph states |12l | as well as 
studying the distinguishability of states by separable op- 
erations [l3|- An important property of the tensor rank 
is that it cannot increase under SLOCC: 



Proposition 1. [l^ Lf 

rkim. 



SLOCC 



IV') then rk{\(j))) 



> 



Through Proposition [U the monotonic nature of the 
tensor rank makes studying it physically worthwhile. Un- 
fortunately, determining the rank of an arbitrary state 
is a very difficult problem [Tsj l which is ultimately why 
there is no simple convertibility test applicable to all tri- 
partite transformations. However, in some special cases 
it is possible to calculate the tensor rank or at least de- 
termine some useful bounds. In this Letter, we establish 
our main results described above by examining the ranks 
of certain tripartite states. We prove the following where 
each statement is in one-to-one correspondence with the 
main results stated earlier. 
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Lemma 1. 

(a') \(t>) G Ha®Hb®Hc can he SLOCC converted from 
state ^I]j=i \^)A\i)B\i)c if and only if rk{\<f>)) < 
N. 

(b') rk{\W)®'^) < 8. 
(c') rfc(|$3)) = 7. 

(d') rfc(|$'^)®") is the multiplicative complexity for mul- 
tiplying two 2" X 2" matrices. 

These results immediately indicate that when many 
copies of a state are considered, the tensor rank does not 
necessarily scale proportionately. As a result, impossible 
SLOCC transformations between individual sates may be 
possible when bulk quantities are considered. 

Extending Lemma 1 to prove Theorem 1 is straight- 
forward. It follows from item (a') that, given a 
tripartite tensor and a number fc, deciding if 
rk(|(/)}) < k can be reduced to the question of whether 

J2i=i\^)A\i')B\i)c ^^-^^'^ 10). The former problem is 
shown to be NP-hard by Hastad [l^, thus the latter is 
also NP-hard (item (a)). 

Results (b), (c), and (d) follow directly from applying 
(a') to (b'), (c'), and (d') respectively. The 17 to 6 con- 
version ratio of (c') is important because 6 copies of 1$'^) 
is a total of 18 EPR pairs. Thus, the stochastic EPR 
distillation rate from multiple copies of \GHZ) is greater 
than 1. In fact (d) shows that this rate can be further 
improved as the upper bound for to is lowered. However, 
the distillation is specific in that the EPR pairs must be 
shared among all three parties. Indeed, if the EPR pairs 
are held by just two parties, r-fe(|$)'*") = 2" so the EPR 
distillation rate from n copies of \GHZ) equals 1. The 
related problem of EPR distillation from the W state has 
recently been studied in [l^. There, the authors show 
that for a single W state, the probability of extracting an 
EPR state via LOCC is not only higher if one does not 
specify which two parties share the state, but it can also 
be made arbitrarily close to one. 

From (d) and the lower bound on /^(2", 2"), it follows 
that 2n copies of GHZ cannot be converted into n copies 
of 1$'^) with a nonzero probability. This result is stronger 
than the one derived in [l3| where the authors demon- 
strate the impossibility of \GHZ)'®'^" |$3^«.n ^^^^^ 
deterministic LOCC. What is most interesting is that 
the authors prove the impossibility strictly through en- 
tropy arguments. Here, we obtain the same conclusion 
using tools of algebraic complexity theory. On the sur- 
face these two lines of attack appear to be unrelated, but 
the similarity in both results suggests that the two may 
be deeply connected. 

Now we turn to prove Lemma 1. We will work with 
unnormalized states below since any overall factor does 
not affect the tensor rank. For any G Ha ^ Hb <E) 



He, let pAB denote Alice and Bob's subsystem obtained 
by taking the partial trace Trc{\(l>){4>\)- As pab is a 
positive operator, it has a spectral decomposition pab = 
SfeLi-Pfe|'0fe)(''/'fe| where < < 1. The vector span 
of {IV'fe) : 1 < < rn} is called the support of pab and 
denoted by supp{pAB)- To proceed, we need the following 
simple equivalent characterization of a tripartite state's 
tensor rank. 

Lemma 2. Suppose € Ha O Hb He- The tensor 
rank of 10) equals the minimum number of product states 
in Ha ® Hb whose linear span contains the support of 
PAB = Trc{\(t)){(f>\)- 

Proof. Let k denote rfc(|0)). Suppose that the span of r 
product states {|aj)|/3j) ■ ^ < j < r} contain supp{pAB)- 
Let 10) = X]i;=i \i)AB\i)c be a Schmidt decomposition of 
10). Each \i) AB belongs to supp{pAB) and thus \i)AB = 
X]j=i ^i,j\'^j)\P3)- Regrouping the \i)c according to the 
r product states gives r > k. On the other hand, consider 
a "minimal" decomposition |0) = X^iLi \(J'i)\bi)\ci)- Then 
PAB = Z]fj=i Wi)\f}i){cj\ci){aj\{bj\ and hence k>r. □ 

Using Lemma [2l the general procedure for determin- 
ing tensor rank is now straightforward. Write |0) = 
J2iLi \i)AB\i)c where the {\i)c : 1 < 1 < m} are or- 
thonormal and then determine the minimum number of 
product states needed to contain the {|«)ab : 1 ^ * £ "m}- 
This question can be rephrased in another way by map- 
ping each \i)AB to a bilinear form from the ring of 
indeterminates (^[{oj}, {6j}] where each aj (bj) is in a 
one-to-one correspondence with a basis vector from Ha 
(Hh). Product states in Ha <E) Hi, correspond to a prod- 
uct of linear forms from C[{aj}] x C[{&j}] which we refer 
to as a non-scalar multiplication. Thus, we obtain the 
following fact: 

Fact. The minimum number of product states that con- 
tain the {\i)AB : 1 < * < '^i}, and hence the tensor rank 
of 10), is the same number of non-scalar multiplications 
Mk ~ (Sj=i '^k,jO,j) X {J2'jLi Pk.jbj) needed to calculate 
the {fi '■ 1 < i < m}. 

We now use the technique outlined above to study the 
tensor rank of certain tripartite states. 

Proof of Lemma 1. (a'): For X^iLi l*)^N)s|*)c'! the 
support of PAB is spanned by N product states. Thus by 
Prop. 1 and Lemma 2, a necessary condition for the given 
transformation is rfc(|0)) < N. Now suppose that |0) = 
J2i=i \a^)\b^)\ci) where k < N. Since {\i)A : 1 < i < N} 
is an orthonormal set, we can define the linear operator 

{la ) 1 < z < A; 
' ' ~ ~ . Similarly, operators B 
0, k < i < N 

and C can be constructed. As noted in the existence 
of such operators is sufficient for an SLOCC protocol 
since |0) will be obtained when Alice performs the local 
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measurement {jj^ 
Bob and Charlie 



I A — ii^ip At A} and similarly for 
Note that (unnormalized) \GHZ)®''^ 



can be expressed as X]i=i W^I*)^!*)^- 
verify by direct computation that \W)'^ 



|ll)A|00)B + |10)A|01)i 



f |01)A|10)i3 

-(|10)^|00)s 
-(|01)a|00}s 



(b'): One can 
expands as: 

|00)a|11)b)|00)c 



|00)a|10)b) 
|00)a|01)b) 
|00)a|00)b) 



|01)c 

|10)c 

|ll)c. 
(2) 



The structure of becomes more manageable when 



working with its corresponding bilinears fi since they can 
be succinctly expressed through the matrix multiplica- 
tion 
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/oi 
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bio 
\biij 



We make use of the following identity 
means a entry: 



where a 
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-aoo/ 



(3) 



Note that rank one matrices require only one non-scalar 

^a^[hi + 62)^ 



multiplication: 


( ai 


ai\ 




Kai 


a J 




any n x n diagonal matrix requires 


Ai \ 




{bi\ 




(Xibi\ 


A2 




b2 




A262 


A3 




bs 




A3 63 


V 




\b,J 




^A464/ 



aj(&i + 62) 



while 



Hence, a total of 



eight non-scalar multiplications is sufficient to compute 
each fi. These multiplications correspond to product 
states that contain supp{Trc{\W) {W\^'^)). By Lemma 
2 then, rk{\W)'^'^) < 8. In fact, expansion ^ gives the 
eight product states that contain supp{Trc{\W) {W\'^'^)) 
enabling us to rewrite Alice and Bob's vector attached to 
\i)c ■ i G {00,01, 10, 11} in (21) as a combination of these 
eight states. To our knowledge, this is the first observed 
non-additivity of the Schmidt measure for pure states, 
(c'): Up to a local unitary transformation on Alice's part, 
the corresponding bilinear forms of |$^) match the set of 
polynomials obtained when multiplying two 2x2 matri- 
ces: 



aoo aoi\ / 600 ^01 
aio an / I 610 611 



foo foi 
fio fll 



(4) 



An algorithm for obtaining the fi using only seven mul- 
tiplications was discovered by Strassen [l9| and later 
proven to be optimal by Winograd [2^ . These seven non- 
scalar multiplications correspond to a minimum number 



of product states containing supp{Trc{\^^) {^^\)) and so 
rfcd^'^)) = 7. As an eight term expansion for liy)"*^ 
can easily be obtained from expansion ([3]), it is straight- 
forward to find a seven term expansion of 1$"^) from 
Strassen's algorithm given in [l3|. Since the explicit ex- 
pressions are not of primary interest here, we omit the 
calculations, (d'): By taking multiple tensor products of 
the matrices in ([4]), we see that for ?? copies of |$'^); the 
corresponding polynomials are represented by 2" x 2" ma- 
trix multiplication. Hence, rfc(|<I>'^)®") is the complexity 
of this operation. □ 

In conclusion, we have found that no easy test ex- 
ists for determining whether two general tripartite states 
are probabilistically convertible. The difficulty arises be- 
cause any general solution involves a tripartite tensor 
rank computation. As a result, one must consider tri- 
partite transformations on a case-by-case basis. In this 
letter we have done this for special states in which the 
tensor rank has already been studied or can be calcu- 
lated with mild effort. Performing this analysis led to an 
improved GHZ state to W state SLOCC transformation 
rate as well as the first demonstration of obtaining EPR 
pairs from GHZ states at a rate greater than one with a 
nonzero probability. 

The connection between tensor rank and entanglement 
transformation is perhaps most beautifully exemplified 
by the equivalence of matrix multiplication complexity 
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and the optimization of stochastic EPR distillation from 
many copies of GHZ. This relationship opens many av- 
enues of further research as the techniques of algebraic 
complexity theory might teach us more about the nature 
and limitations of SLOCC transformations. Conversely, 
constructing explicit SLOCC entanglement transforma- 
tions using results in quantum information may be useful 
to obtain bounds for the multiplicative complexity of a 
particular set of bilinear forms. 
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